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ABSTRACT 

The  transition  to  three-dimensionality  in  a  plane  mixing  layer  is 
studied  through  exploration  of  the  flow  field  and  analysis  of  the  power 
spectrum  of  longitudinal  velocity  fluctuations. 

Two  different  phenomena  are  identified.  First  the  appearance  of  3-D 
small  scale  fluctuations  is  marked  by  the  transition  of  the  slope  of  the 
inertial  subrange  towards  the  value  of  -5/3.  Second  we  find  in  one  case 
that  the  large  scale  itself  shows  an  instability,  with  wavelengths  of  the 
order  of  the  vorticity  thickness  of  the  layer,  and  that  seems  to  be  expla¬ 
ined  by  a  collective  "corrugation"  ot  the  vortex  cores. This  second  pheno¬ 
mena  is  studied  in  some  detail. 
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1.  INTRODUCTION 

Ic  is  generally  agreed  that  three-dimensionality  is  an 
essential  property  of  fully  developed  turbulence.  Even  if 
2-D  turbulence  is  possible  it  is  only  under  three-dimension¬ 
al  conditions  that  a  fully  developed  energy  cascade  is 
established  and  leads  to  the  production  of  enough  small 
scale  motion  to  account  for  the  rapid  mixing  of  concentra¬ 
tion  and  momentum  gradients  that  is  characteristic  of 
turbulent  flees. 

The  plane  mixing  layer  is  a  well  studied  flow  in  which 
the  initial  conditions  are  two-dimensional  and  the  flow 
remains  two-dimensional  for  a  long  time  after  separation.  In 
fact  there  is  evidence  that  the  large  scales  of  the  flow  may 
stay  two-dimensional  essentially  forever  [1,2].  Since  it  is 
obvious  that  the  small  scale  notion  in  the  nixing  layer 
should,  and  indeed  does,  eventually  become  three-dimension¬ 
al,  this  flow  is  a  good  test  case  to  study  the  instabilities 
that  are  responsible  for  the  breakdown  of  two-dimensional 
motion. 

This  transition  has  been  studied  previously  by  various 
investigators  from  several  points  of  view.  Bradshaw  [3] 
measured  the  evolution  of  the  turbulence  level  in  the  mixing 
layer  and  ita  approach  to  its  equilibrium  value.  Perhaps 
significantly  this  approach  was  found  to  depend  on  whether 
the  initial  boundary  layer  was  turbulent  or  laminar  and 
thus,  possibly,  on  the  amount  of  initial  three-  dimensional¬ 
ity  in  the  flow.  The  evolution  of  the  turbulence  level  has 
been  measured  more  recently  by  other  investigators  for 
different  initial  conditions  [4,5]  but  it  is  still  not  clear 
whether  the  establishment  of  three-dimensionality  coincides 
with  the  attainment  of  the  equilibrium  value  of  the  fluctua¬ 
tion  intensity. 

More  recently  Konrad  [6]  and  Breidenthal  [71  have 
identified  a  transition  in  the  amount  of  molecular  mixing 
present  in  a  two-species  shear  layer.  Since  enhanced  mixing 
must  come  from  an  increase  of  small  scale  structure  this 
transition  is  probably  associated  with  the  onset  of 
three-d imensinality  [8]. 

Me  have  used  still  another  property,  the  evolution  of 
the  alope  of  the  inertial  subrange  of  the  longitudinal 
velocity  power  spectrum  towards  its  3-D  value  of  -5/3,  as  an 
indicator  of  three-dimensionality.  The  transition  is 
clearly  visible  in  the  three  cases  that  we  have  studied  and 
is  discussed  in  section  5. 

Another  interesting  and  related  problem  is  the  possible 
appearance  of  three-dimensionality  in  the  large  structures, 
as  opposed  fu  the  small  scales.  Chandrsuaa  et.  al.  [ 9 1 
reported  helical  distorsion  of  the  vortices  in  a  plane 
mixing  layer  and  similar  effects  have  been  described  in  the 
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initial  region  of  circular  jets  [10].  Kourad  has  reported 
the  appearance  of  longitudinal  streaks  in  his  concentration 
pictures  and  Breidenthal  has  associated  these  streaks  with 
sinusoidal  corrugations  of  the  vortex  cores.  In  one  of  our 
flows,  the  only  one  checked  for  detailed  apanwise  structure, 
we  have  found  pernanent  longitudinal  jets  that  nay  corres¬ 
pond  to  the  sane  effect  and  seen  to  be  associated  to, 
although  they  precede,  the  transition  to  snail  scale  three- 
dinensionality .  These  results  are  described  in  section  6. 

Sons  of  our  neasurenents  were  carried  out  with  hot  film 
anenoneters  in  water.  Since  that  technique  is  not  too  well 
developed  a  significant  aaount  of  tine  was  spent  chocking 
and  inproving  the  accuracy  of  the  method,  and  the  results 
are  briefly  discussed  in  section  3. 

Due  to  the  tine  limitations  involved  in  this  study  many 
of  the  experiments  presented  are  exploratory  and  many  of  the 
conclusions  tentative.  Further  work  is  being  carried  out  to 
confirm  and  extend  then. 
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2.  FLOW  FACILITIES  AND  EQUIPMENT 

The  experiments  have  been  performed  in  cwo  independent 
facilities!  a  water  tunnel  and  2-D  air  half-jet. 

Water  Tunnel 


The  experimental  apparatus  is  basically  a  blow-down 
tunnel  in  which  a  plane  mixing  layer  is  established  between 
two  streams  of  water  coming  from  two  large  independent 
reservoirs  (figure  1).  Both  tanks  are  pressurised  with  air 
at  1.3  kg/ca?  to  keep  the  change  in  velocity,  due  to  the 
drop  in  water  level  throughout  the  shot,  below  a  maximum  of 
.52 . 

The  cwo  streams  go  through  a  conventional  sequence  of 
control  valves,  flow-meters.  settling  chamber,  straightnars 
and  screens  and  a  4:1  contraction  before  they  are  brought 
into  contact  to  form  a  plane  mixing  layer  at  the  exit  ot  two 
7.5  x  13  cm  nozzles  in  the  test  section,  from  where  they  are 
discarded  through  a  siphon. 

The  test  section  is  IS  x  15  cm  across  and  100  cm  lung. 
In  the  lateral  walls  there  are  three  pairs  of  registers 
which,  together  with  a  sliding  probe  holder,  allow  manual 
positioning  of  the  probe  at  any  downscream  coordinate  up  to 
4u  cm  from  the  edge  of  Che  splitter  plate.  Cross-stream 
probe  motion  ia  achieved  by  means  ot  a  traversing  gear  which 
incorporates  a  stepping  motor,  with  a  .2  mm  step,  and  is 
controlled  by  a  HP  2100S  minicomputer  which  automatizes  the 
water  experiments. 

Although  Ll.t  apparatus  can  provide  a  flow  with  tempera¬ 
tures  as  well  as  velocity  differences,  the  experiments 
reported  in  this  work  have  been  conducted  in  a  homogeneous 
shear  layer. 

A  steady  flow  is  established  in  less  than  a  second  at 
uominal  velocities  of  19  and  37  cm/s.  Under  these  conditions 
the  momentum  thickness  of  both  boundary  layera,  ac  cue  end 
of  the  splitter  plate,  are  .7  and  .3  mro  respectively.  A 
waiting  time  of  about  15  s  is  needed  before  the  water 
initially  in  the  settling  chamber  is  swept  completely  and 
the  temperature  is  uniform  enough  to  allow  the  use  of  hot 
film  sensors.  With  this  in  mind,  maximum  useful  running  time 
is  about  50  3. 

We  have  run  performance  tests  on  the  apparatus  and  the 
free  stream  turbulence  level  was  found  to  be  about  .2%.  The 
free  stream  velocities  were  found  repeatable  for  100  differ¬ 
ent  shots  and  constant  throughout  a  given  one  to  within  .52. 
The  shear  across  the  central  part  of  the  free  stream  high 
velocity  side  was  less  than  .22  and  of  the  order  of  12  in 
the  low  velocity  side. 

Even  if  the  tanks  are  thermally  Insulated  it  is  unnvoid- 
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able  to  have  soma  ctnptracurt  difference  between  them.  Thia 
difference  ia  monitored  by  meana  of  thermistors  and  can  be 
kept  below  .1  *C  which  repreaenta  a  velocity  error  of  about 
.21. 

Air  Half  Jet 

The  air  flow  facility  was  used  during  Che  last  part  of 
the  axper lisente .  A  blower-tunnel  was  utilized  to  supply  a 
tvo-d iaensionsi  nixing  layer  in  air,  A  blowing  fan  powered 
by  a  3.7  Kw,  2870  rpm  aotor  supplied  the  air  flow  through  a 
diffuser  to  a  planus  chamber ,  The  plenum  chaaber  contained 
foas,  screens  and  a  deep  honeyconb  in  order  to  reduce  the 
turbulence  level  at  the  noaxle  exit.  The  air  was  exhausted 
through  a  rectangular  noxale  of  7x13  era  (figure  2.).  The 
contraction  ratio  waa  10:1. 

The  jet  was  allowed  to  six  on  one  of  its  boundaries  with 
the  surrounding  quiescent  air,  the  ether  three  being  solid 
to  prevent  Che  formation  of  mixing  layers. 

A  traversing  gear  provides  .hree  degrees  of  freedom  with 
a  precision  of  .05  mm  in  the  x - y  plane  and  .2  mm  in  the 
spanvise  direction. 

The  maximum  valocity  could  be  aa  high  as  50  m/s  although 
all  the  experiments  in  air  have  been  performed  at  25  or  1 6 , ^ 
m/s,  A  more  detailed  description  of  this  facility  is  givo 
in  section  6. 

Anemometry  and  Data  Acquisition  System 

Velocity  measurements  were  carried  out  in  both  mix  >»j 
layers  with  a  DISA  55D01  anemomoter  in  constant  tempera' u.e 
mode  without  linearization.  The  output  of  the  anemometer  was 
fed  into  a  10  bit  HP5610A  ana  log -d  ig i t a  1  converter  and  the 
digitized  data  were  stored  by  the  HP  21Q0S  minicomputer  in 
magnetic  tapa.  Data  rate  waa  governed  bv  mi  independent 
programmable  crystal  clock  which  provided  a  s..  ’’le  frequency 
standard.  The  taps  thus  produced  was  then  processed  off-line 
in  a  large  IBM  370/158  computer. 

In  the  water  tunnel  facility  the  ex^riif  *nts  were 
controlled  by  the  HP  2100S  that  opened  .me1  closed  the 
valves,  moved  the  probe  and  monitored  suh  ,i ic  i  i  ' .-»  * 1  data  as 
tank  temperature,  flowmeter  output  and  p'o’oe  pootion. 
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3.  MEASUREMENT  TECHNIQUES 
Ho  ■  Film  Anemomt  r  y  in  Water 

The  velocity  measurements  in  water  were  carried  out 
uaing  a  hot  film  quartz  coated  single  eenaor  DISA  55A81 
wedge  probe.  Probe  overheat  was  kept  at  .05,  which  gave  us 
enough  sensitivity  while  minimizing  the  risk  of  bubble 
formation  and  deposition  of  dirt  on  the  probe  surface  [111. 
Even  at  thia  low  overheat,  and  even  if  we  filtered  our  water 
through  a  3  u«  filter  to  minimize  the  number  of  nucleation 
points,  bubbles  were  present  (see  figure  3)  and  about  10Z  of 
the  runs  had  to  be  discarded  because  of  this  problem. 

Surface  contamination  produced  a  day  to  day  drift  in  the 
thermal  characteristics  of  Che  film  which  caused  a  change  in 
Che  calibration  curve.  To  study  this  effect  we  calibrated 
our  probe  in  a  small,  low  spaed  recirculating  tunnel  vith 
very  low  turbulence  level,  in  which  the  velocity  was 
measured  simultaneously  vith  the  probe  and  with  a  calibrated 
flowmeter  placed  in  the  recirculating  duct.  The  upper  pert 
of  figure  4  shows  different  calibrations,  carried  out  on 
different  days  and  showing  clearly  the  drift.  To  correct 
this  effect  we  tested  several  heat  transmission  models  based 
on  the  expression 

Nu  -  Nu  -  i»NuU  (  3 .  1  ) 

r 


where  Nu  is  the  Nusaelt  number  as  actually  measurad  by  the 
probe,  Nu,.  i~  the  corrected  Nusaelt  number  that  would  be 
measured  by  an  ideal  uncontaminated  film  and  fl  is  a  parame¬ 
ter  related  to  the  ratio  of  thermal  conductances  of  the 
coating  and  the  fluid  that  must  be  found  empirically  for 
each  calibration.  A  detailed  account  of  this  method  can  be 
found  in  F121.  There  we  find  that  the  best  collapse  of  the 
data  happens  for  n-2  (aae  figure  5)  in  agreement  with  the 
model  proposed  by  Morrow  and  Kline  in  fl3‘.  Applying  this 
correction  to  the  calibration  data,  all  points  fall  close  to 
the  "ideal"  curve  that  ia  shown  in  the  lower  part  of  figure 
4.  The  standard  deviation  with  respect  to  this  curve 
decreases  from  8Z  before  correction  to  L6S, 

The  application  of  this  correction  model  to  experimental 
situations  presents  additional  problems,  since  at  least  the 
real  velocity  at  one  point  has  to  be  known  to  compute  i  . 
Diaphragm  flowmeters  were  introduced  in  the  supply  lines  to 
monitor  the  volumetric  flows  but,  since  Che  test  section 
walls  were  fixed  and  we  expected  some  pressure  gradient  to 
develop  on  the  acreatawice  direction,  these  flows  gave  os 
only  the  velocity  at  the  exit  of  the  nozzles.  In  fact  the 
pressure  gradient  had  to  be  computed  us  part  of  the  probe 
calibration  procedure. 
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Lat  ua  assume  chat  one  data  point  ia  measured  In  «ach 
stream  giving  Nusselta  Nuj  and  Nu_>  and  Chat  the  velocities 
at  Cha  axle  of  the  nozzles  arc  known  to  be  Ojp  and  .  The 
velocities  in  the  sections  under  consideration  can  be 
related  Co  the  pressure  through  Bernouilli's  aquation  which 
after  linearization  gives, 

”i  '  ’’iE  "  »PA*"iE  0.2, 

If  the  calibration  formula  for  tha  proba  ia 

l»  -  f(Nu  )  -  f(Nu  -sNu2),  (3.31 

linaarizing  again  and  substituting  in  (3.2)  we  get  from  the 
two  streams  s  system  of  linear  equations  permitting  ua  to 
compute  6  and  tbs  presaura  as  functions  of  the  measured 
Nuaaalt  numbers  and  valocitiea  at  the  nozzle.  Figure  6  shows 
tha  pressure  gradient  computed  in  this  way  in  five  different 
co-ord inates  along  Cha  teat  section,  confirming  tha  consis¬ 
tency  of  tha  above  procedure  and  tha  applicability  of  the 
correction  method. 

Taking  into  account  the  scatter  and  errors  associated 
with  the  calibration  of  tha  probe,  the  flowmeters  and  the 
msaaureaents  of  water  temperature  and  pressure  gradient  we 
estimate  tha  accuracy  of  our  velocity  measurements  to  be 
approximately  31, 

Hot  Wire  Measurements  in  Air 

In  the  mixing  layer  in  air  classical  hot  wire  nnemometry 
techniques  were  applied.  Velocity  luctustions  were  measured 
with  a  DISA  55F31  hot  wire  probe,  ’  urn  in  diameter  am)  1 .25 
ram  in  length.  An  overheat  ratio  of  0.8  was  used  in  all 
cases. 

Data  processing 

Two  different  types  of  data  were  taken  during  the  water 
tunnel  experiments:  a)  profile  measurements  at  four  down¬ 
stream  sections  and  b)  spectral  measurements  at  fixed 
positions  within  the  shear  layer. 

Ten  to  tvelve  points  were  taken  in  each  traverse.  The 
data  collection  period  per  point  was  3  seconds  which  corres¬ 
ponds  to  approximately  10  large  structures  passing  by  at  the 
most  downstream  section,  x*37  cm. 

The  spectral  measurements  were  taken  leaving  the  probe 
standing  still  at  fixed  positions  within  the  layer  and,  in 
order  to  resolve  very  low  frequencies,  the  data  collection 
period  at  these  points  was  the  whole  duration  of  the  run, 
i.e.  40  to  50  seconds. 
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In  the  air  experiments  the  measurements  taken  were 
essentially  the  same,  a* though  the  exploration  was  extended 
"o  th-ee  dimensions,  and  tue  data  collection  peifiod  was 
changed  to  correspond  to  the  different  frequencies  involved. 

Mean  and  rms  velocities  in  water  were  computed  by  using 
the  calibration  curve  on  the  digital  data,  and  taking  the 
average  and  si &ndar 3  deviation  of  the  resulting  velocity 
records.  Our  previous  knowledge  of  the  free  stream  veloci¬ 
ties^  Che  pressure  gradient  and  the  hcc  film  drift  was  used 
in  : his  reduction. 

The  spectra,  both  in  air  and  v/ater,  were  also  computed 
in  tnis  way.  The  mean  velocity  wa's  removed  for  each  run,  and 
an  averaged  FPT  technique  [14],  with  a  Blackraann-Harris  data 
window  [15],  wac.  -..feed  to  obtain  the  power  spectra.  Finally 
all  the  spectra  c< rrespod ing  to  a  given  station  were  aver¬ 
aged  together.  Since  the  le.igth  of  a  long  run  is  40  to  50  s 
tid  4  to  6  runs  were  made  on  every  single  point  the  minimum 
number  of  large  structures  presented  in  each  averaged 
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4.  RESULTS 


Flow  description:  water 


The  top  part  of  figure  7  shows  nond ioenoional  profiles 
of  mean  velocity  which  fit  fairly  well  to  an  error  function 
curve.  In  spite  of  the  slight  favourable  pressure  gradient 
mentioned  above  the  growth  is  linear  and  similarity  is 
achieved.  The  vorticity  thickness  of  the  layer  was  measured 
and  the  spreading  rate  found  to  be  .0793,  with  a  virtual 
origin  located  26  ou  upstream  of  the  splitter  plate  edge. 
This  spreading  rate  is  perhaps  a  little  high  but  within  the 
range  of  values  given  in  previous  studies  [16],  [17].  The 
dividing  streamline  is  na-«0065  where  the  nondimensional 
mean  velocity  is  (U-U2) / (U iwU2 ) “ • 58 . 

The  bottom  part  of  the  figure  presents  the  rms  value  of 
longitudinal  velocity  fluctuations  at  two  sections.  In  the 
most  upstream  one  some  influence  of  the  initial  conditions 
can  still  be  detected  and  the  general  turbulence  level  is 
higher  than  in  the  other  one.  The  Reynolds  numbers  for  both 
sections',  xaU/v,  are  34000  and  63000  which  are  below  the 
limits  given  in  [3]  and  [10]  as  needed  to  attain  similarity 
of  the  fluctuation  profiles.  In  fact  figure  9  shows  the 
maximum  value  of  fluctuations  measured  as  a  function  of 
downstream  coordinate  and  it  is  clear  that  it  has  not  yet 
attained  its  equilibrium  value.  This  asymptotic  value  can  be 
estimated  to  be  between  .18  and  .19  which  is  again  a  little 
high  but  reasonable.  The  facts  that  both  the  growth  rate  and 
the  turbulence  level  are  slightly  high  might  nor  be  inde¬ 
pendent  [18], 

Some  time  lapse  pictures  of  the  shear  layer  were  taken 
by  injecting  dye  at  the  splitter  plate  and  these  pictures 
were  used  to  measure  the  visual  spreading  rate.  The  result¬ 
ing  value  <$ / (x-xq ) ■ . 1 25  is  consistent  with  that  given  in 
[19],  even  if  it  is  not  clear  whether  dye  visualization 
should  give  the  same  spreading  rates  as  the  shadowgraphs. 


Flow  d escr iption : 


air 


A  study  of  the  two-dimensional  air  half  jet  was  under¬ 
taken  with  the  splitter  plate  of  type  I  of  figure  2  ,  which 
is  essentially  a  thin  aluminium  sheet  and  a  jet  velocity  of 
25  m/s.  The  mixing  layer  thus  formed  has  been  analyzed 
through  mean  and  fluctuation  velocity  profiles  and  spectra 
taken  at  several  positions. 

The  top  of  figure  8  presents  nondimensional  profiles  of 
mean  velocity,  which  fit  fairly  well  to  the  expected  error 
function  curve.  From  these  profiles  we  see  that  similarity 
has  been  attained  at  least  since  x»80  ram.  The  vorticity 
thickness  of  the  layer  was  measured  and  the  spreading  rate 
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was  found  to  be  .195,  with  a  virtual  origin  located  9  am 
upstream  of  the  splitter  plate  edge.  This  spreading  rate  is 
half  w ay  between  those  reported  by  Wygnansky  and  Fiedler 
[20]  and  Liepmannand  Laufer  [21],  The  dividing  streamline  is 
110*703  and  the  corresponding  non-dimensional  mean  velocity 
is  (U-U2 ) / (Ui -U2 )■ . 57 .  This  result  is  the  name  as  the  one 
found  in  [ 21 ] . 

The  bottom  part  of  figure  8  shows  the  rms  value  of 
longitudinal  velocity  fluctuations  at  three  sections.  The 
most  upstream  profile  indicates  that  the  turbulent  intensity  * 
has  not  achieved  its  equilibrium  value,  and  this  is 
confirmed  in  figure  14.  The  asymptotic  value  of  the  maximum 
u*/AU  is  .17  which  again  agrees  fairly  veil  with  those 
reported  in  [20]  and  [21]. 

A  third  experiment  using  the  same  facility  and  a  veloci— . 
ty  of  16.4  m/s  is  described  in  section  6, 


5.  SPECTRAL  RESULTS 

In  order  tc  stud/  the  development  and  interaction  of  the 
large  coherent  structures  and  their  possible  relation  to  the 
origitl  and  evolution  of  three-dimensionality  ve  have  carried 
out  a  detailed  spectral  analysis  of  the  velocity  fluctua¬ 
tions  in  our  three  mixing  layers:  the  two-stream  water 
facility  and  the  air  half-jet  at  two  different  velocities. 

All  spectra  were  taken  at  the  mid-span  plane.  Those  in 
water  were  measured  either  at  n “y/ (x-*o )“0  or  n».Q43.  The 
last  position  permitted  a  reliable  detection  of  the  large 
structure  peaks,  not  yet  swamped  by  the  strong  low-frequency 
noise  that  seems  to  be  associated  to  the  center  of  the 
nixing  layer,  while  still  shoving  a  high-frequency  inertial 
subrange  characteristic  of  fully  turbulant  flow.  A  typical 
evolution  of  the  power  spectra  across  the  mixing  layer  is 
shown  in  figure  10. 

For  the  high-dpeed  air  jet  (2Sm/a)  all  spectra  were 
taken  at  n-0  while  those  for  the  low  speed  one  (16.4  m/s) 
fall  in  two  series,  none  of  them  along  a  similarity  ray.  One 
series  was  measured  at  U/Um^.S  and  was  used  to  study  the 
inertial  subranges  while  the  other,  at  U/U^^.OS  and  there¬ 
fore  in  the  outer  edge  of  the  layer,  was  used  for  the 
detection  of  the  large  structures.  Since  this  work  is 
related  to  mixing  layer  transition  and  most  of  the  data  are 
far  from  the  asymptotic  similarity  range,  this  later 
approach  was  found  more  satisfactory  then  keying  the 
measurements  to  geometric  variables  that  are  relevant  mostly 
for  the  downstream,  well  developed,  part  of  the  layer. 

The  frequency  behaviour  of  the  spectra  can  be  divided 
broadly  into  three  ranges:  low  frequencies,  large  structures 
and  inertial  range.  Only  the  last  two  are  studied  separately 
below  since  the  low  frequency  one,  though  interesting  C22l, 
does  not  seem  to  be  influenced  by  three-dimensionality. 


The  spectrum  of  the  large  eddies 

The  part  of  the  spectrum  corresponding  to  the  large 
structures  is  the  most  interesting  one  but  also  the  most 
complex  as  can  be  seen  in  figure  II,  which  shows  the  down¬ 
stream  evolution  of  the  water  spectra. 

One  of  the  most  noticeable  aspects  is  the  existance  of 
peaks  which  can  be  followed  across  several  stations  with  a 
very  definite  constant  frequency,  while  the  low  frequency 
ones  gain  importance  as  we  move  downstream.  The  evolution  of 
the  position  of  these  peaks  with  downstream  coordinate  is 
given  in  the  bottom  part  of  figures  12-a,b  and  c,  for  the 
three  layera.  For  the  purpose  of  these  figures  a  peak  was 
included  if  its  maximum  power  was  at  least  102  above  the 
spectral  continuum  at  that  frequency. 
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An  important  phenomenon  is  that  the  ratio  between 
successive  spectral  lines  is  mostly  two,  although  different 
factors,  like  1,5,  are  also  present.  These  values  agree  with 
those  given  by  Usman  and  Jimenez  [23]  and  Jimenez  et  al. 
[22  ]  and  can  be  explained  by  means  of  the  amalgamation 
process.  The  doubling  is  easy  to  understand  through  the 
coallescence  of  two  equal  eddies.  The  factor  l.S  ia  also 
explained  assuming  the  amalgamation  between  structures  of 
two  successive  generations,  which  implies  chat  some  eddies 
of  the  last  one  have  escaped  without  pairing  and  hence 
coexist  and  eventually  amalgamate  with  vortices  whose  scales 
are  twice  as  large. 

The  same  variables  are  shown  nond imensionalized  in 
figure  13.  The  frequencies  are  made  non-dimensional  with  the 
convection  velocity,  UaU1-fU2/2,  and  with' a  character  1st ic 
lenght  which  is  supposed  to  be  proportional  to  the  wavel¬ 
ength  of  the  initial  instability  of  the  layer.  For  the  air 
jets  this  length  is  the  momentum  thickness,  62  *  the 
boundary  layer  [3,  19];  for  the  water  layer  we  have  used  an 
equivalent  thickness,  related  to  the  most  unstable  wave 
length  of  the  wake  profile  of  the  splitter  plate  (see 
Appendix  A).  The  longitudinal  coordinate  ia  made  non-dimen¬ 
sional  with  the  length  U<52/AU,  proportional  to  the  length 
needed  for  a  spatially  amplified  disturbance  to  develop 
C24]. 

In  these  dimensionless  variables  the  frequencies  of  most 
spectral  peaks  agree  for  the  three  mixing  layers,  especially 
near  the  splitter  plate,  indicating  that  instability  theory 
dominates  the  flow  for  at  least  a  few  pairings.  Those 
spectral  peaks  which  appear  only  in  one  of  the  layers  can 
usually  be  traced  to  forcing  frequencies,  e.g.  blower 
frequency,  which  are  not  intrinsic  to  the  flow.  Moreover,  in 
the  downstream  region  in  which  some  degree  of  self-preserva¬ 
tion  has  been  achieved,  all  layers  follow  a  unique  frequency 
law  which  can  be  written  as 

A  =  .3  AU  x/U  (5.1) 

and  agrees  very  well  with  Che  data  reported  by  other  inves¬ 
tigators  C25,26]. 

Some  coupling  between  the  coherent  structures  and  the 
instabilities  of  the  initial  boundary  layer  in  the  air 
experiments  is  apparent  in  that  frequency  peaks  that  will 
later  be  dominant  in  the  mixing  layer  are  visible  very 
early  in  the  downstrem  development  (fig.  12-b,c).  The 
boundary  layers  in  these  cases  are  marginally  unstable 
(U  650  and  520)  ,  and  some  of  the  frequencies  induced 
either  by  the  blower  fan  or  possibly  by  coupling  with 
downstream  structures  are  amplified  and  seem  to  survive  in 
the  mixing  layer.  Other  frequencies  which  are  far  from  the 
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instability  region  in  the  boundary  layer  do  not  survive  and 
are  lost  quickly,  even  if  their  absolute  amplitude  is  larger 
(see  also  figure  16  and  discussion  in  section  6).  In  the 
water  mixing  layer  the  boundary  layers  at  both  sides  of  the 
plate  are  stable  (U61/v<  450)  and  the  flow  is  "cleaner"  than 
in  the  air  jets. 

The  inertial  subrange 

Jill  the  spectra  in  figure  10,  at  different  transverse 
positions,  show  an  inertial  subrange  with  a  slope  of  -5/3; 
those  in  figure  11,  taken  at  various  downstream  stations, 
also  show  the  same  slope,  except  the  one  taken  at  x«150  no 
which  is  much  steeper,  close  to  -3.  Since  -5/3  is  associated 
to  the  three-dimensional  energy  cascade  the  study  of  this 
subrange  can  be  used  to  detect  the  outset  of  threed imensior.- 
slity . 

The  upper  parts  of  figures  12-a,b  and  c  present  the 
slopes  of  the  inertial  subranges  in  all  those  stations  in 
which  it  could  be  measured  with  some  confidence.  It  can  be 
seen  that  the  behaviour  of  the  experiments  in  air  and  voter 
is  quite  different.  The  water  flow  (figure  12a)  presents  a 
relatively  sharp  evolution  from  a  slope  of  -3  to  -5/3,  The 
value  of  -3  is  associated  to  the  vorticity  cascade  which  is 
normally  assumed  to  control  two-dimensional  turbulence  [271, 
and  this  suggests  that  the  observed  transition  separates  a 
situation  in  which  velocity  fluctuations  are  predominantly 
two-dimensional,  even  in  the  small  scales,  from  one  in  which 
3-1)  fluctuations  become  important.  The  change  of  slope 
occurs  in  a  region  in  which  several  scales  of  coherent 
structures  co-exist  (see  bottom  of  figure  12a)  implying  that 
Che  flow  is  more  disordered  there  than  in  other  zones. 

"he  mixing  layer  of  the  air  experiments  do  not  show  any 
evidence  of  a  plateau  near  the  slope  of  -3  and  their 
appr jach  to  -5/3  is  more  gradual.  Since  the  turbulence  level 
i.  the  free  stream  of  this  facility  was  somewhat  higher 
(1%),  and  the  boundary  layers  sore  unstable,  the  origin  of 
tur eu-Uimensionality  could  be  governed  by  different  causes 
M  in  the  water  flow.  In  fact  the  situation  is  somewhat 
reminiscent  of  the  behaviour  of  the  maximum  turbulence  level 
in  layers  in  which  the  initial  conditions  are  laminar  or 
turbulent  [3,  4]  , 

Figure  14  shows  the  slopes  and  the  maximum  turbulence 
levels  for  the  three  layers  in  terms  of  the  same  non-diaen- 
. 'tonal  coordinate  used  for  the  sp^ctr&l  peaks.  Two  facts  can 
ba  seen  in  this  figure.  Firac  the  transitions  in  water 
happens  sooner  and  is  indeed  sharper  Than  the  one  in  air. 
This  is  again  consistent  vjf  the  idea  that  initial  tu~Su- 
lance  tends  to  retard  tran  . ;  v.ion.  Second  the  approach  to 
three-dimensionality  is  not  completed  until  the  velocity 


fluctuations  begin  to  decrease  to  their  asymptotic  level, 
confirming  that  the  tvo  phenomena  are  not  unrelated. 

The  high  frequency  end  of  the  inertial  range  is  marked 
by  a  corner  in  the  spectrum  beyond  vhich  the  slope  becomes 
much  steeper.  The  position  of  the  corner  is  given  in  the 
bottom  parts  of  figures  12a, b,c  and  is  very  near  in  all  cases 
to  the  inverse  of  the  sensor  length,  no  that  the  new  range 
is  probably  explained  by  the  limitations  of  the  sensing 
equipment . 

Finally  a  word  must  be  said  concerning  the  dimensionless 
variables  used.  Considering  that  only  three  independent 
experiments  are  available  it  is  not  really  possible  to 
decide  in  most  esses  which  is  the  appropriate  normalization 
and  the  variables  have  been  chosen  with  some  particular 
theoretical  model  in  mind.  The  growth  of  the  large  struc¬ 
tures  is  the  best  understood  process  and  consequently  the 
non-dinenjionalizatiou  works  best  in  that  case  (figure  13). 
The  cause  of  the  three-dimensionality  is  less  clear  and  we 
have  done  little  elae  than  trying  several  normalizing 
factors.  Using  ad  aheissa  the  ratio  x/<S2  gave  slightly 
better  results  than  the  one  used  here,  but  was  difficult  to 
justify  theoretically  for  the  tv  l-stream  case.  The  obvious 
parameter  , AU  S^/v,  increased  the  differences  between  the 
three  experiments  to  almost  an  order  oi  magnitude. 
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6.  SPANWISE  FLOW  STRUCTURE 

An  important  question  is  whether  the  transition  to 
three-dimensionality  is  itself  two-dimensional.  An  obvious 
possibility,  for  example,  is  that  perturbations  are  injected 
from  the  end  plates  that  limit  the  shear  layer  in  span,  and 
this  mechanism  has  been  shown  to  be  important  in  some  cases 
of  transition  in  the  two-dimensional  wake  behind  a  cylinder 
[28],  If  this  be  so  in  our  case  it  should  be  expected  that 
the  three-dimensionality  will  appear  further  upstream  as  we 
approach  Che  end  plates  and,  to  test  this  hypothesis,  we 
have  explored  the  flow  in  a  shear  layer  along  the  z  coordi¬ 
nate  as  well  as  in  the  x?y  plane.  The  somewhat  surprising 
results  are  presented  in  this  section. 


Experimental  arrangement . 


The  measurements  were  carried  out  in  the  air  half  jet 
adjusted  to  s  nominal  velocity  of  16.4  m/s.  Two  slightly 
different  splitter  plate  configurations  were  used  which  are 
shown  in  outline  in  figure  2.  Cronologically ,  configuration 
1  was  used  first  and  later  taken  appart  and  substituted  by 
configuration  II.  The  boundary  layer  in  both  cases  was 
laminar  with  a  momentun  thickness  of  about  .2  mm  at  the  exit 
section,  corresponding  to  an  effective  1-snght  of  flat  plate 
of  90  mm  in  case  I,  and  slightly  less  in  case  II.  This 
thickness  was  both  computed  from  the  area  distribution  of 
the  exit  nozzle  and  measured  directly  by  exploration  of  the 
flow  with  a  hot  wire  just  downstream  of  the  plate  edge.  When 
severe  spanwise  nonuniformities  in  the  flow  were  discovered 
in  shear  layer  I  the  two-dimensionality  of  the  edge  was 
checked  and  found  to  have  irregularities  with  an  amplitude 
of  nearly  .5  mm  and  a  characteristic  wavelength  of  50  mm. 
Also  a  slight  wobble  (.2  mm)  was  found  in  the  y  position  of 
the  probe  as  it  was  traversed  in  the  z  direction.  As  a 
consequence,  and  even  if  it  was  doubtful  whether  such  smooth 
irregularities  could  cause  the  observed  effects,  the  split¬ 
ter  plate  was  substituted  by  a  carefully  machined  square 
prism  (II)  whose  two-dimensionality  is  thought  to  be  better 
than  .1  mm  in  the  130  mra  span.  The  traversing  mechanism  was 
also  rebuilt  to  minimize  the  wobble  and  any  remaining 
deviations  were  calibrated  against  the  edge  of  the  splitter 
plate.  This  calibration  procedure  vas  repeated  several  times 
and  the  resulting  positional  accuracy  of  the  probe  with 
respect  to  the  plate  edge  was  estimated  to  be  about  .05 
mm(rms).  The  x  and  y  movement  of  Che  probe  was  done  by  a 
pair  of  crossed  lathe  beds  whose  positional  resolution  is 
also  . ? 5  mm,  which  probably  accounts  for  much  of  the  uncer¬ 
tainty  quoted  above.  The  data  reported  in  this  section  refer 
to  layer  II  unless  explicity  stated  otherwise. 
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A  general  sketch  of  the  nozzle*  geometry,  together  with 
the  coordinate  system  used  and  the  the  approximate  position 
of  the  mixing  region  is  also  ahovn  in  figure  2. 

All  measurements  ware  taken  with  a  single  hot  wire  1.25 
bb  long  and  5  pm  in  diameter,  whose  holder,  with  a  diameter 
of  3  mm,  was  held  at  an  angle  of  45  degrees  out  of  the  plane 
of  the  flow  to  minimize  interference.  The  wire  itself  was 
aligned  parallel  to  the  splitter  plate  edge. 

The  spanvise  uniformity  of  the  free  stream  was  checked 
and  found  to  be  better  than  .51.  The  residual  deviations 
were  smooth  and  showed  no  sign  of  oscillatory  behaviour.  The 
free  stream  was  2X  faster  near  th*  splitter  plate  than  at 
the  center  of  the  jet  in  a  way  consistent  with  the  expected 
flow  at  the  sxit  of  a  concentration.  The  level  of  longitudi¬ 
nal  velocity  perturbations  <in  the  stream  was  comparatively 
high  (.81)  and  concentrated  in  a  few  well  defined  frequen¬ 
cies,  mainly  harmonical  of  tha  rotation  frequency  of  the 
blower  fan  (48  Hz),  and  as  a  consequence,  the  boundary  layer 
contained  a  rather  high  level  of  perturbations  (4Z) .  Power 
spectre  for  the  free  stream  and  boundary  layer  longitudinal 
velocities  are  given  in  figure  16  and  show  a  clear  corres¬ 
pondence  between  Che  important  peaks  in  both  eases.  The 
frequencies  of  these  peaks  are  also  plotted  on  the  Toll- 
mien-Lin  stability  diagram  for  the  instability  limit  for 
two-dimensional  perturbations  T291  and  those  peaks  at  the 
approximate  frequencies  show  signs  of  slighe  amplification. 
Although  no  measurements  were  taken  of  the  mechanical 
vibrations  of  the  splitter  plate  it  was  clear  y  feeling 
manually  tha  aide  walls  of  the  tunnel  that  pai  t  of  the 
perturbations  to  the  boundary  layer  were  transmitted  that 
way,  and  by  taking  a  series  of  steps  to  change  the  rigidity 
and  mass  of  the  tunnel  walls  the  level  of  perturbations  was 
reduced  from  about  10Z  to  the  final  4Z. 

The  boundary  layer,  measured  just  downstream  of  the 
splitter  plate,  was  not  very  two-dimensional.  Four  profiles 
are  given  in  figure  17  in  which  the  transverse  coordinates 
have  been  normalized  with  the  local  measured  momentum 
thickness.  The  average  velocitie  are  seen  to  collapse 
satisfactorily  and  to  follow  qm.  .  closely  tha  Blaaius 
profile  for  the  laminar  layer,  but  the  turbulence  levels 
show  a  large  scatter.  Figure  18  shows  the  spanwise  variation 
of  the  momemtum  thickness  and  Che  maximum  turbulence  level 
across  half  of  the  total  span.  The  momentum  thickness  shows 
variations  of  almost  20Z  and  the  non-uniformity  of  the 
turbulence  level  is  even  higher.  It  is  not  clear  which  is 
the  cause  of  this  nonuniformity  although  it  is  presumably 
connected  either  with  small  geometric  irregularities  in  the 
nozzle  walls  or  wl'.h  changes  in  the  amplitude  of  the  mechan¬ 
ical  vibrations.  The  fact  chat  the  layer  is  just  transition¬ 
al  probably  amplifies  the  effect  of  any  such  perturbations. 
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The  raea sur er.ien C 8  vara  done  over  a  period  of  several  days  and 
enough  points  were  re-checked  to  make  sure  that  the  effect 
was  repeatable  and  that  the  position  of  the  nonuniformities 
remained  constant  with  time,  but  no  detailed  comparison  was 
done  between  the  two  splitter  plate  configurations  and  it  is 
not  known  which  was  the  effect  of  rebuilding  the  plate  on 
the  boundary  layer. 

Results 


Most  of  the  results  presented  in  this  section  refer  to 
average  properties  integrated  analogically  over  periods  of 
about  one  minute  and  reduced  to  velocity  without  making  any 
allowance  for  effects  due  to  hot  wire  nonlinearity.  An 
exception  is  the  data  used  in  power  spectra,  which  were 
digitized  directly  and  treated  in  the  way  described  earlier 
in  the  report.  Each  spectrum  derives  from  about  400  a  of 
measurements. 

Using  splitter  plate  I,  we  took  data  at  two  different 
spanwise  locations,  mid— span  (z-0)  and  z«35  mm.  Average 
longitudinal  velocities  and  turbulence  levels  corresponding 
to  the  mid-span  position  are  shown  in  figure  19  where  the 
transverse  coordinate  i3  made  nond imens iona 1  with  the  local 
computed  vorticity  thickness  and  normalized  to  the  local 
central  streamline,  defined  as  the  point  where  U-.5  Uro.  It 
is  seen  that  the  profiles  become  seif  similar  beyond  x*60 
mm,  and  that  this  is  true  even  for  the  turbulence  levels, 
although  with  a  larger  scatter. 

A  less  reassuring  view  is  presented  in  figure  20  where 
the  growth  of  the  vorticity  thickness  and  the  evolution  of 
the  maximum  turbulence  levels  are  given  for  three  cases,  the 
two  spanwise  locations  in  layer  I  and  the  mid-span  position 
of  layer  II.  It  is  specially  interesting  to  note  that  the 
growth  of  the  thickness  of  layer  l  is  not  the  same,  even  at 
the  most  downstream  coordinate,  for  the  two  spanwise  loca¬ 
tions,  and  that  the  data  for  the  other  configuration  do  not 
agree  either  with  any  of  the  two  curves,  even  if  the  three 
asymptotic  growth  rates  are  within  the  range  of  values  found 
by  other  investigators. 

As  a  consequence  of  this  result  the  nozzle  configuration 
was  rebuilt  and  a  detailed  spanwise  exploration  of  layer  II 
was  undertaken.  Figure  21  shows  the  turbulence  levels 
measured  on  the  plane  ym— . 082  (x-3 . 45  mm),  and  it  is  clear 
that,  also  in  this  case,  there  are  spanwise  nonuniformities 
amount  iug  to  nearly  13%  of  the  average  turbulence  level  and 
that  these  nonun  i  f ormit ie s  organize  themselves  broadlv  in 
long  longitudinal  streaks. 

The  next  three  figures  show  maps  of  average  velocities 
and  turbulence  levels  in  three  roughly  orthogonal  sections 
of  the  shear  layer.  Figure  22  is  a  conventional  x-y  map 
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taken  at  mid-span  and  shows  ti.«  expected  velocity  distribu¬ 
tions.  In  particular  the  peak  that  appears  near  x*30  mm  in 
the  turbulence  level  seems  to  be  characteristic  of  flows 
with  a  laminar  initial  boundary  layer  r 3 1  . 

The  other  two  figures  are  spanwise  sections  and  are 
constructed  in  a  slightly  different  way.  For  each  x-y 
position  a  representative  value  of  velocity  and  turbulence 
level  ia  defined  by  averaging  over  the  spanwise  coordinate, 
x,  and  the  figures  show  the  difference  between  actual 
measured  properties  and  those  locally  constructed  averages; 
they  represent  the  effect  of  spanwise  nonunif oraities 
independently  of  the  general  background  levels  of  the 
velocity  and  turbulence  fields. 

The  x-z  section  in  figure  23  shows  clearly  the  stream- 
vise  turbulence  streaks.  The  corresponding  average  velocity 
map  shows  that  these  streaks  are  associated  to  long  stream- 
vise  jets  of  higher  than  normal  velocity.  The  jets  appear  to 
attain  their  full  strength  near  x«40  mm  and  are  involved 
later  in  a  series  of  interactions  that  keep  increasing  their 
separation  progressively,  although  discontinuoualy,  without 
appreciably  decreasing  their  intensity. 

The  transversal  structure  of  these  streaks  is  shown  in 
the  y-z  section  in  figure  24  where  the  jeta  are  shown  to 
have  transversal  as  veil  as  longitudinal  coherence.  Perhaps 
the  most  significant  feature  in  chi*  figure  is  that  the 
velocity  field  is  roughly  simmetrical  about  the  central 
streamline,  while  the  turbulence  level  is  approximately 
antysimmetrical.  This  is  in  contradiction  with  the  explana¬ 
tion  of  the  nonuniformities  as  streamwise  vertices  of  the 
type  postulated  by  Benney  [30]  and  known  to  be  important  in 
the  transition  of  the  laminar  boundary  layer,  and  agrees 
more  with  a  model  involving  atreamwise  jets,  or  cross-strea;; 
vortices.  This  will  be  discussed  later  in  more  detail. 

In  understanding  these  jets  an  obvious  possibility  is 
that  they  are  somehow  inherited  from  some  irregularity  in 
the  initial  conditions.  We  saw  before  that  the  free  stream 
in  the  neighbourhood  of  the  splitter  plate  was  exarainated 
and  shown  to  be  uniform  within  the  accuracy  of  the  hot  wire 
equipment,  but  chat  the  spanwise  nonunif ormities  found  in 
the  boundary  layer  just  after  separation  were  important 
enough  to  be  considered  as  likely  candidates  Co  be  the  cause 
of  the  stronger  features  found  further  downstream. 

before  this  matter  is  examined  a  word  must  be  said 
regarding  the  accuracy  of  the  measurements.  It  was  noted 
above  that  Che  position  of  the  probe  was  known  only  to 
within  .05  mra  and,  since  near  the  splitter  plate  there  are 
strong  transverse  velocity  gradients,  this  uncertainty 
induces  errors  in  the  measurement  that  can  be  important  in 
the  upstream  part  of  the  layer.  These  errors  are  large 
enough  that  data  on  velocicy  variations  upstream  of  x-I5  mm 
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are  not  really  significant  and  could  be  consistent  with  a 
uniform  average  velocity  while  turbulence  data,  being  less 
affected  by  strong  gradients,  are  qualitatively,  although 
not  quantitatively,  reliable  all  the  way  to  the  splitter 
plate . 

Having  in  mind  this  warning,  the  x-z  map  in  figure  23 
suggests  that  there  ia  no  direct  continuity  between  the 
features  in  boundary  layer  and  those  downstream.  lv.  particu¬ 
lar  the  region  between  x“5  and  20  no  seems  to  be  one  of 
reorganization  of  the  spanwise  structure,  while  the  fact 
that  the  separation  of  thr  jets  changes  downstream  and  that 
there  is  a  rapid  increase  ia  thair  intensity  nser  x»40  mo 
suggests  chat  some  kind  of  local  instability  of  the  mixing 
layer  is  responsible  for  them. 

If  we  think  of  the  large  coherent  eddies  of  the  mixing 
layer  as  representative  of  the  dominant  longitudinal  mode  of 
layer  instability,  the  spanviBe  nonunif'ormit  lee  can  be 
considered  as  a  secondary  wave  superimposed  on  that  main 
one,  and  by  comparing  the  power  spectra  of  the  flow  at 
varioua  z  positions  it  sfould  be  possible  to  study  the 
properties  of  that  wave. 

For  a  n articular  x-v  position  we  define  a  reference 
curbulenco  level  «u  -  and  power  spectrum  P  by  averaging 
over  the  span.  At  a  given  %  the  deviations  from  these  values 
can  be  ascribed  to  the  secondary'  spanwise  wave  and,  to  a 
firat  approx imat ion ,  we  could  expect  that  the  deviation  in 
the  spectra  are  proportional  to  those  in  a^.  If  this  is  true 
the  normalized  excess  power  spectrum 


p  ( f ,  z )  -  p  ( n  n 

p  (  t' ,  Z  )  *  - - — T~  (  b  .  1  ) 

P  (  f )  e  •  (  z  )  -  o  ‘ 

should  be  independent  of  z  and  represent  the  frequency 
structure  of  the  secondary  wave,  Figure  25  shows  an  estimate 
of  this  quantity  for  throe  different  x-y  points,  obtained  by 
evaluating  equation  (6,1)  at  several  spanwise  positions  and 
averaging  over  z, 

p  (  f  )  =  *  P  (  f .  7.  >  “•  2  (  b  .  2  ) 

•> 

Since  the  integral  of  the  power  spectrum  should  equal  o' 
it  is  easy  to  see  that  p  should,  in  Che  average,  be  equal  to 
1,  with  those  frequencies  where  p>l  representing  the  spec¬ 
tral  p< aks  associated  with  the  secondary  wave.  The  main 
oeaks  evident  in  the  spectra  in  figure  25  are  represented  in 
the  lower  part  of  the  figure  superimposed  on  the  x-fraquenev 
diagram  that  gives  the  evolution  of  the  large  structure  in 
the  layer  (see  also  figure  12-c).  If  the  jets  were  just 
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"fossil"  effects  deriving  from  activity  in  the  boundary 
layer  one  would  oxpect  their  spectral  behaviour  to  remain 
constant  downstream,  with  peaks  Associated  either  to  the 
apnnwlse  instabilities,  of  the  transitional  boundary  layer 
[31]  or  vith  the  forcing  frequency  from  the  blower  fAn. 
Instead.  the  dominant  frequency  shifts  downstream  and 
remains  locked  approximately  with  the  frequencies  associated 
to  the  large  eddies  in  the  mixing  region.  This  gives  further 
support  to  the  view  chat  the  screaks  are  caused  by  an 
intrinsic  local  instability. 

A  final  piece  of  evidence,  although  fragmentary  and 
difficult  to  interpret,  is  important  enough  chat  it  should 
be  presented  here.  Before  the  apparatus  for  layer  I  was 
taken  appart  and  the  exit  noezle  rebuilt,  some  exploratory 
measurements  were  taken  of  the  spanwise  distribution  of 
velocity  and  turbulence.  Since  it  way  expected  that  the 
reworking  of  the  splitter  plate,  with  its  complete  change  in 
the  position  and  nature  of  all  random  irregularities,  would 
change  also  the  nonuniformities  in  the  flow,  no  measurements 
were  taken  for  layer  II  at  the  same  points  as  in  I. 
Nevertheless  the  measurements  taken  in  I  are  represented  in 
figure  ? 6  together  with  their  closest  equivalent  in  II. 
While  the  most  upstream  data  do  not  show  any  similarity  in 
both  cases,  due  may  be  as  much  to  the  inaccuracy  of  the 
measurements  in  that  station  as  to  anything  else,  the 
situation  is  different  in  the  downstream  sections  where  the 
s imil ar it  it s  are  high  enough  Chat  it  is  difficult  to  escape 
the  conclusion  that  the  two  flows  have  basically  the  same 
spanwise  structure. 

Since  there  is  a  very  low  probability  that  the  spanwi.se 
structure  of  the  initial  boundary  layer  is  even  roughly 
similar  in  both  cases,  the  similarity  of  the  final  flow 
seems  to  imply  that  the  initial  conditions  are  not  important 
for  the  development  of  spanwise  nonunif o rra i t ios .  Whether 
this  persi9tance  of  the  structure  after  the  reworking  of  the 
splitter  plate  is  due  to  some  intrinsic  property  of  the 
layer  or  to  an  irregularity  of  the  free  stream,  not  the 
boundary  layer,  small  enough  to  be  undetectable  by  our 
equipment  but  still  important  for  the  flow  ia  not  clear  and 
must  await  further  investigation. 

Theoretical  considerat ions 

We  have  shown  that  under  some  conditions  a  mixing  region 
that  should  be  expected  to  be  two-dimensional  contains  a 
definite  three-dimensional  structure  consisting  of  high 
turbulence  streaks  associated  to  high  velocity  longitudinal 
jets.  The  position  of  these  streaks  remains  fairly  constant 
from  one  day  to  another  and  even  shows  some  degree  of 
permanence  after  severe  changes  in  the  initial  conditions. 
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Similar  three-dimensional  affects  have  b«an  observed  in 
other  f  lovs  ,  notably  in  the  tranaitional  boundary  layer  f 3 1 J 
and  iu  the  vaka  behind  a  two-dimensional  body  [28  1,  In  the 
mixing  layer  Konrad  [6 J  observed  longitudinal  streaks  in  the 
concentration  interface  while  Breidenthal  [7]  repeated  the 
observation  under  different  conditione  and  showed  that  the 
streaks  arise  f.roa  “corrugations'*  in  the  vortex  core*  <  f  the 
large  structure. 

Benney  [30]  studied  theoretically  the  nonlinear  deve.l«.  ■■•» 
ntnt  of  the  instability  of  a  laminar  mixing  layer  anrl 
concluded  that  a  three-dimensional  vavn  cim  Intnrucr  with 
ths  dominant  three-dimensional  inode  Co  produce  permanent 
three  dimensional  variations  of  average  velocity  anaoelated 
to  longitudinal  vorticas  alignad  with  the  x-axia.  Tb*.t  work 
was  aimed  at  axplaining  tho  transition  in  rhc  ;.««inar 
boundary  layer,  where  longitudinal  vorricea  are  knew-  to 
apperr,  but  unfortunately  does  not  explain  the  tranavereal 
structure  of  the  streaks  observed  in  our  case. 

In  fact  a  longitudinal  vortex  transfers  momentum  accrosa 
the  mixing  layer  producing  an  average  antys iramet r ic  distri¬ 
bution  of  velocity  excess,  in  contradiction  with  tha  results 
presented  in  figure  24. 

Tha  structure  shown  in  thst  figure  suggest  that  tha 
vorticas  responsible  for  the  streaks  are  aligned  in  the 
transverse,  y,  direction,  while  the  spectral  results  show 
that  these  vertices  srs  somehow  in  phase  with  the  iarge 
structure.  When  the  two  results  are  put  together  we  arrive 
at  a  picture  in  which  che  vortex  cores  of  the  primary 
two-dimensional  instability  are  deformed  or  corrugated  in 
accordance  with  Breidenthal ' s  pictures,  and  the  corrugation 
has  at  least  a  component  in  the  y-r  plane.  At  Che  points 
where  the  corrugation  is  concave  towards  the  high  velocity 
side  the  velocity  is  higher  while  che  center  of  che  core, 
and  therefore  the  region  of  higher  turbulence,  is  displaced 
towards  the  low  speed  side  of  the  layer  (see  figure  27). 
This  picture  is  in  qualitative  agreement  with  che  data  in 
figure  24  but,  if  these  features  are  to  be  visible  in  the 
average,  all  the  cores  have  to  deform  “in  step". 

Saffman  T32]  has  suggested  that  the  instability  respon¬ 
sible  for  che  corrugation  of  tho  cores  is  the  same  one 
appearing  in  vortex  rings.  Briefly,  the  criterion  for 
instability  is  that  a  core  deforms  in  such  a  way  that  its 
self-induced  rotational  speed  io  zero  [33],  in  which  case 
the  strain  field  induced  by  other  cores  in  the  layer  is  able 
to  pull  it  appart  and  feed  the  growth  of  che  deformation. 
This  growth  starts  at  an  angle  of  45  degrees  to  the  flow 
direction  and  has  therefore  the  transverse  component  implied 
by  our  observations. 

The  lowest  mode  for  which  this  happens  is  a  corrugation 
of  the  core  in  which  both  radial  and  spanwise  vorticity 
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distributions  are  changed.  The  moat,  unstable  wavelength 
depends  on  Che  vorticity  distribution  assumed  initially 
within  the  cores  but  it  falls  between  2  and  2.5  .  tines  the 
vortex  radius  [34,  35],  and,  if  we  estimate  this  radius  as 
about  half  the  vorticity  thickness  of  the  layer  we  arrive  at 

A  *  1  to  1.25  6U  (6.3) 

While  this  is  only  a  rough  estimate,  the  wavelengths 
measured  in  the  layer  and  shown  in  figure  28  are  indeed  of 
this  order  of  magnitude,  specially  in  the  downstream  part  of 
the  layer  where  the  streak  system  is  well  developed.  Thus, 
this  mechanism  might  be  able  to  explain  the  instability  of  a 
single  vortex  core,  but  the  question  of  whether  it  is 
possible  to  lock  all  the  vortices  in  the  layer  into  a 
coherent  instability  chat  will  produce  the  observed  jets  is 
still  open  and  needs  further  investigation. 


7.  CONCLUSIONS 


There  seems  to  be  at  least  two  different  phenomena 
connected  with  the  three-dimensional  breakdown  of  the  plane 
mixing  layer.  First,  at  least  in  some  cases,  there  is  a 
macroscopic  instability  of  the  large  structure  which  produc¬ 
es  a  spanwise  system  of  longitudinal  jets  whose  amplitude 
builds  up  very  quickly  and  then  remains  approximately 
constant  while  their  separation  adjust  itself  to  a  value 
that  is  somewhat  higher  that  the  vorticity  thickness  of  the 
layer.  This  phenomenon  may  be  due  to  a  collective  instabili¬ 
ty  of  the  vortex  cores  of  the  same  type  as  the  one  found  in 
vortex  rings. 

The  second  phenomenon  is  the  development  of  a  3-D 
inertial  subrange  in  the  spectrum  of  the  small  scales  of  the 
velocity  fluctuations.  In  the  only  case  in  which  we  have 
documented  both  phenomena  the  second  one  doea  not  develop 
fully  until  the  macroscopic  instability  has  not  been  growing 
for  some  time  and  has  actually  arrived  to  what  seems  to  be 
its  final  configuration. 

Moreover  there  is  some  indication  that  there  are  two 
qualitatively  different  ways  for  the  establishment  of  the 
3-D  inertial  ranges,  exemplified  here  by  the  water  and  air 
experiments.  In  the  first  one  the  transition  is  abrupt  and 
there  is  an  initial  part  in  which  the  inertial  range  has  a 
slope  characteristic  of  the  2-D  enstrophy  cascade.  In  the 
second  one  the  transition  is  more  gradual  and  there  is  no 
evidence  of  a  2-D  regime.  It  is  not  clear,  at  present,  which 
is  the  condition  responsible  for  the  differences  between  the 
two  cases. 

An  interacting  coincidence  is  that  in  the  three  cases 
studied  here,  the  trsnsition  to  three-dimensionality  in  the 
small  scales  happens  near  the  second  or  third  vortex  amalga¬ 
mation  (see  figures  12-a,b  and  c),  despite  differences  of 
one  order  of  magnitude  between  the  Reynolds  numbers 
involved.  It  may  be  that  the  amalgamations  act  as  trigger 
for  the  transition,  but  the  small  number  of  cases  makes  it 
impossible  again  to  derive  definitive  conclusions. 

We  believe  the  work  presented  in  this  report  should  be 
seen  aa  a  preliminary  exploration  of  the  effects  connected 
with  the  3-D  transition.  Much  more  research  is  needed  in  all 
cases.  In  particular  it  is  important  to  check  whether  the 
spanwise  macroscopic  breakdown  can  be  reproduced  in  differ¬ 
ent  circumstances  and  how  it  depends  on  initial  conditions 
and  Reynolds  number.  An  interesting  question,  in  view  of  the 
importance  of  3-D  structure  for  a  vigorous  mixing,  is 
whether  this  instability  can  be  enhanced  and  used  to  shorten 
the  distance  needed  for  transition  in  practical  flows. 

Secondly  the  influence  of  vortex  amalgamations  on 
transition  should  be  checked.  Again  more  experiments  with 
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different  initial  conditions  are  needed  for  this,  but  the 
tools  developed  in  this  study  using  the  power  spectra  of 
velocity  fluctuations  should  be  enough  to  detect  both  the 
average  position  of  the  pairings  and  the  position  of  three- 
diaensional  breakdown. 

We  want  to  acknowledge  th<*  support  of  part  of  this 
investigation  by  the  Spanish  Cot  :  sioa  Asesora  de  Investiga¬ 
tion  as  well  as  the  help  of  the  UAM-IBM  Scientific  Center  in 
Madrid  which  gave  us  access  to  its  computer  center  for  the 
digital  processing  of  the  data  in  this  report. 
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APPENDIX  A 

Stability  characteristic a  of  the  splitter  plate  wake . 

the  wavelength  of  the  initial  instability  in  a  nixing  layer 
is  strongly  influenced  by  the  wake  of  the  splitter  plate  in 
those  cases  in  which  both  free  stream  velocities  are  differ¬ 
ent  froa  zero.  The  downstream  distance  needed  for  viscosity 
to  fill  the  dip  in  the  profile  produced  by  the  initial 
boundry  layers  is  of  the  order  of  the  equivalent  flat  plate 
length  of  those  layers  while  the  amplification  distance  for 
the  unstable  disturbances  is  proportional  to  the  layer 
thickness  and  is  usually  auch  smaller.  Moreover,  for  the 
coamonly  found  Reynolds  numbers  the  wake  flow  is  unstable 
froa  the  moment  it  leaves  the  splitter  plate  and  the  nonli¬ 
near  regiae  of  the  instability  is  achieved  before  any 
appreciable  evolution  of  the  profile  has  taken  place. 

It  seeas  likely  therefore  that  the  wavelength  of  the 
initial  instability  is  fixed  by  the  initial  wake  profile  and 
we  present  here  an  analysis  of  a  simple  straight  line 
approximation  to  this  profile.  Approximation  of  this  type 
have  been  shown  to  be  useful  in  other  cases  [36,  37]  and 

give  a  quantitative  basis  for  comparison  between  different 
profiles.  In  particular  the  variation  of  wavelengths  between 
different  velocity  ratios  should  be  given  approximatly  by 
this  model. 

Tha  profile  we  uae  is  shown  in  figure  29.  If  we  intro¬ 
duce  a  perturbation  stream  function 

VCx,  y,  t)-Y(x,y)+4>(y)  expfia(x-ct)  )  ,  (A.l) 

and  assume  infinite  Reynolds  number  and  parallel  flow,  :he 
perturbation  amplitude  satisfies  the  equation  [38], 

(U-c)(V-a2<t)-U"$=0.  (A. 2) 

In  the  segments  in  which  U  is  linear  the  solution  of  (A. 2) 
is  given  by 

<t^exp  (i  ay),  (A  .  3) 

while  at  the  corners,  where  U"  is  infinite,  the  stream 
function  remains  continuous  but  the  derivatives  satisfy  a 
jump  condition 

4>'  -$>'  =  $  (U  '  ~  U  '  )  /  (U-c  )  (A. 4) 

t"  “  *T“  “ 

When  we  add  the  conditions  that  t>  remains  bounded  at  infini¬ 
ty  and  assuming  that  the  real  part  of  a  is  positive,  the 
solution  of  (A.l)  can  be  divided  in  the  segments  shown  in 


figure  29.  The  jump  conditions  give  relations  between  the 
various  coefficients,  and  the  final  result  is  a  cubic 
eigenvalue  equation  to  be  satisfied  by  c. 

In  the  case  of  the  nixing  layer  the  equivalent  flat 
plate  length  for  both  boundary  layers  is  usually  the  sane, 
and  the  thickness  are  related  by 

.  h2/h!  -  (U1/U2)1/2  ( A . 5 ) 

The  solutions  to  the  eigenvalue  equation  with  this  condition 
is  given  in  figure  30,  where  the  stability  boundary  ia  given 
in  terns  of  the  wave  number  and  the  velocity  ratio,  together 
with  isolines  of  the  spatial  anplif ication  rate  ah^C£/cr. 
The  properties  of  the  nost  amplified  wave  are  also  given  in 
figure  31  as  functions  of  the  velocity  ratio. 

For  our  water  mixing  layer  U2/Ui«.51  and  the  most 
unstable  wave  number  ia  1.25  times  chat  of  the  half  jet 
with  the  same  boundary  layer  thickness  at  the  high  speed 
side.  Therefore  the  wavelength  and  the  "equivalent"  initial 
thickness  can  be  taken  as  .8  times  the  real  high  speed  side 
value,  which  is  the  thickness  used  in  section  5. 

Some  of  the  ideas  presented  in  this  appendix  where 
suggested  by  Prof.  A.  Linar;. 
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Fig.  1. 


PRESSURIZATION 


-  Schematic  view  of  water  installation. 


Variation  of  rms  error  after  correction  of  calibration 
drift  as  a  function  of  the  exponent  in  equation  O.H. 


(> .  -  Pressure  gradient  in  test  section  of  water  tunnel 
computed  as  part  of  the  drift  correction  procedure 
for  the  hot  f i In. 
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Fig.  7.-  Average  velocity  (top)  and  fluctuation  profiles 
for  the  water  mixing  layer;  51  . 
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Average  velocitv  C cop)  and  fluctuation  profiles  for  the  air 
haif-jet  (Uoo=25  m/s). 


Fig-  12a,  b,  c.-  Evolution  of  the  frequency  of  the  spectral  peaks 

corresponding  to  large  scales  (o) ,  and  the  frequency 
of  the  "nose”  marking  the  end  of  the  inertial  range 
(A).  Also  slope  of  the  inertial  range,  a)  Water  layer; 
b)  Air  (25  m/s);  c)  Air  (16.4  m/s). 
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Fig.  13.-  Dimensionless  representation  of  large  scale  frequencies  versus 
downstream  coordinate.  O  i  water;  X  ,  Air  (25  m/s);  £  ,  Air 
(16.4  m/s) . 
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Fig.  15. ~  Definition  of  the  inertial  and  orobe- limited  ranges  in  the 
power  spectrum,  showing  the  "nose"  separating  them. 


17.-  Nondimens  tonal  profiles  of  boundary  laver  at  several 
spanvise  locations.  O  »  z  -  -35  mm;  X  »  z*-25;V  , 
z  ■  15;  A  .  z  =■  -5 .  Solid  line  is  Bias  ius  profile. 
Laver  IT,  x  •  I  mm. 


18.-  Spanwise  variation  or  maximum  turbulence  level  (— — ) 

and  momentum  thickness  l - 1  tor  boundary  laver. 

I  .aver  II.  x  »  1  mm . 


Fig.  19.-  Nondimensicnal  longitudinal  ve.locicv  and  trubui enc 
level  profiles.  For  trite  air  layer,  o  .  x  =  40  mrii; 

X  *  x  »  60  •  A  ,  x  =  80 ;  7  ,  x  -  1  00’.  Laver  I ,  z  =  0  ; 

Uo»  =  lf>.  4  o/.s. 
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Fig.  26.-  Normalized  deviation  of  power  spectra  from  spanwise 
uniformity;  p  defined  in  text.  Each  curve  is  average 
of  several  z  positions.  x*l,  y**.05  mm:  z  = -40,-30, 

-15,  10,  20,  35  mm;  x*20,  y«*-.9  mm:  z=-35,  -25, 

-15,  5,  10,  25,  40  mm;  x=60,  y  = -6  mm:  z*-40,  -35, 

-25.  -20,  -10,  -5,  5,  15,  25mm.  Spectral  peaks  are 
overlaid  on  frequency  diagram  for  layer,  (see  figure 
12c)  to  show  relation  to  dominant  peaks  in  the  spectr-um. 


Fig.  29.-  Straight  line  approximation  for  the  wake  profile  of 
the  splitter  olate.  Expressions  are  the  perturbation 
amplitudes  for  che  stream  function. 


